


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1978-12 


Pressure on the interface between a 
converging fluid wedge and a fast fluid bottom 


Kawamura, Masami: loannou, loannis 


Monterey, California. Naval Postgraduate School 
http://ndl.handle.net/10945/18522 
Copyright is reserved by the copyright owner 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist sia Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ia) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


’ ’ ¥ 
' d , ’ ; t ee 
4 ait 
. 
| 7 13 } 
’ t ‘ 3 u *¢ 
; i “s 
' 
i) 
o's oy » x, , 
‘ f) 
: 4 Ue 
i) 
4 : y 
a 7 Ne 4 ) 
ia 
4 y* 
i] 
hes 
4 
3 
4 ‘ 
*’ 
' a ‘ 
’ 
HY 
4 1 i 
* 
, 
a 
1 
1 
. ‘ 
3 
} ' 
1 i 
tery 
1 ' 
4 
‘ 
' 
| x 
’ 
4 
-~ 
' 
* 
f 
4 
‘ 
A : 
‘ 
A 
| an 
) 
* 
' 
i 
5 
{ 
ry i 
: cp | 
¥ i 
\ 3 
ri 
’ 
' J 
: i l t 
t : 
i oo 
F 4 
‘ 
1 
, 
Sagi 
' a 
‘ ah 
3 4, 
pes t 
: oe 2 
‘ 
‘ , % ¥*> i 


ae 
' 


wide, ees Je, | ee ee eRe eek 

a ida 4 Mak ¥ CARA NE st vay AN 
5 : Rec ie 
mGie 4 * Hy , ash * 4: % bs 4 






























i PA Re Aa Bae 
4 Zz A RAL RAN ORLA TY bo ANS 
hire ee Bt atten yd AR 
] ie aK an ‘4 AAG rie ay Ky 
ns y ae HA za YY CN ‘ hy ee 
6 ait ta 3 «4% % iaa bial Port's 4,7, 
Vd ae Tg a tai ee SA ae 
th rh wey ie Bh VO any ay" 
ay Uh RY Gat Ree ae aN 
4 eu * ai ae t deol ee 
Aw | ; a] ok ue ai 4 ba aa line 
' Ry j 4 FA} %, iS im wes 
P tre wih a a> : pha » WAR RC ot” eft 
wre Nae vy ‘ath, i oe Tee 
’ \ fraygrecj | 1 AY: a he ura +6 ‘ 
¥a : ‘ R : f } iy va . 7 f 
1 5 th, neh he SBE AAG Y We sy 
" 1% re a “y 5 ae ad 4 org i 4 NSB Te Pa het 
a #8 Le + ; 1u4 ar 4 ‘ a * “ed, ‘y is a‘ Yay Wana ny 
VAT NT MALE age hte on fi acs Aas vine eit ie 
y ca he He a t ¢ , Od ee, +t . 
PGE TEASE Ras YS Se 
: yt ‘ a) Y \ WWNea alt! NAN ans ye nt hs Te aS et 
: ney” ¥ ’ 4 bee « ACO Pie Lanes ay, F ‘S ayy RAT 
a 1 <a “why bas a an! NRG . nel Ve Waite 
a | or ae yl ‘I moe he i ea way aan Se NO wa BAEC 
' 7 o. 7 ee La aw ‘y Sent sy hei aT “yy lay WS i 
ha ny Lae yi vit tr ! h : } 4, CN 1 Rd hel Thy 
ry Mh eZ eg fp HG. LOM uti’ : ° > 
" 7h a 23 94, iit a Me \a\e.$ aay » P FRA en eae bse = 
at: is sa x Ne yt ax ayy ah a eo Bess si 
; ‘. ' J ' . * +h 4 .) * rh v 
; Tere sh sire: > Wee, A IMs 
ay Y ' a4 yl] t Zz Ae 4 7 4 Li e “ ea! : i*)'s te ; 
a . , 4S Ae a ‘ i % ae #18 a: can n ie 
ta * Lt \ af de ‘Pe ply ° Ve, . hee ; k ’ asf: i 5 
“ 1 1 os WN aly \ a ty > ty * etary) sty’ Al GA 
ec 4 a ve > % ga ‘ i & ty ! re F ae 4 Lots on) 
ws Pe : 4 1’) ae Va j a4 ° _ .] 
eo ‘ ye, “A . , Cn Cee | i pay 
P , er ‘ . 4 1 mpi - » thes aay at 
mm \ : , oe is ae a AT egod™ ae ot Se a an oye Nei ay ts 
‘ 2 Yo a. 8 v . * WR em Aaheye e 2 =e ee at ; 
; \ N 4 ‘ rs ASY ; 5° : br ho 7 4 = i Sve a ees = 
: : “ ' . y TF : — aPal thee Fe aay 
, . . Vitae . <5 ey ace Ses EA : in 
7 a x as ‘ ‘ Z ae = F Bain’ 
; \ 4 5 A bd i a bf a rt Wee sete i hee . . ‘x - 
i an Be whiefe % a a 7 ah. ” 
A ‘ ae “ton = ' w = «  & — 
. ; a ‘ ‘ % . a ag * 3 ie ot en “a 
* = = = ble . a 
a4 . % P y * . {s i %e Rw = 
: - i ~ a “i " ' oy ile ee ae ay 3) ro 
‘ Ce bor s »’ a , -- 4 Ros toh es 
i # ,“ . oie ee , ‘ Ls rr. hevem 2 Ths SE Tee, 
is * - Bit Ss ‘ ie L ad ' acts aid 
* Y - - w~ a e. lag | peste a 
e hc Pe et I ein 
io a oy en we X% be os! > aor s 
a 4 = = x —i — i ‘ = = i Any - 
; . ew ? TT. te 2 ie 9 a6 “> . = ne fi 
» = am mt . “-_* « @ a _ SBS. & tte _—~ re 3 
I ten , > = eat « a at athe SUT | Sei eaey eet Al 
i i ; aa ‘= i we. a ak 
A 5- . oor Sn we ata mee 5 -_ . 
a — b/ s ‘ea i Lawn es i 
A % LY a " i ww llae So ae = 
a . 1 - i a ‘ . 
‘ r = j . === F 
4 atg= 4 é Od = & ate! ~~ d 
ty ’ ' : | = _ ie 
i “ = = 
i = _ Pe " i] cS 
1 ~ i ‘ Wael 
. ’ A “ at | . 
a1 P 2 a ae a 
Ge ! ' - 
‘ - ! & hihe - 
| wy 
t J = sa 
: ha ' i ' t a " 
fi 4 ~ : ey Ae 
i i” - é 4 
ay t a a a 
. : a a] J 
¥ im * y me ee 
Jas 3 -_ 
$e \ = e.= 
; : Ta 
? ; he \ = 7 4 + ! Pet 
- i 4 4 r f= 
1a yeas i 
| , c i 
1 ‘ UT 
i =" - Tit = 
rt 
. : , ' J 4.) 45) | > \ ws 
. is ee hee 
e iat hes? 
1 es 2 = 13 tas eet el watts! 
i LU = — ssa 
“sn Fa i. ae ite 
"lee \ i = = | os = 
a a i 
| il a pi f {-e = " 4) i hy 
it j a7 ’ i = a 
t ' L] 
= i | ! oot y 
P j a ‘ Ly 
q L | | mi 
r ; Au: h ! ae | 
i a 
I \ 
e Cie 7 , mt i ' ‘ 
i + 
ce = Wi i 
i) I iv 4 i , r ! poe} é J ' b i = | +¥ 
i i - i? = us 1 - 
i. '; tal . ! i i ** 4. ‘ - i i] i 4 
fg OT { fe a ‘i, ~ J 
“eee Or ] i ti us 
fe ey s ’, J 
‘ : ei r om 
i + & | 
; | | PAL Me pe Ay] 
i 1 i i fn ‘ oes Fe iy Fs - SERIE a - 
Me | 
' ; 
t ri “4 
i fez f 
: Ca ae 
: k a a 
* } *, a ie 
! r F j j ge in ‘ 
; ¥ ay , | y “4 4 me ¢ ey vy 
a 5 a, 1 4 , 4 { ; Ps Py ot - ‘f a & tit 
a ae rte Ue W ! 
1 oi f b 
1 } of t 
4, i Li 
/ i 5 4 fr 4 BY: 93 
hs pe i 4 oF ’ Kobe ‘ f hay Att 
dw’ $ 3 ry MFI G | “A f fe Raa f Tis yap oe 4 Vega 
F , $ for’ Se i ¥ ‘Ss sin yore iy 4 
: EERE CRE ern E Si 
‘ f , vf Ay ry v3 rn # a ¥ ? 
‘ 2 veh Nery Ae neat ee te Ta spe Py sh * 
a) " a bk he x *e ", ai ty ; iad oy pers: oi , 
4 i f pty. i ¢* oe i" ; ay) iy Ln fal Pa wy dy ry 
: pL Nk Meng BP BN aI 
tp er J j i vi j s0, y ahy Me ie ‘ 
t s zy | V1 Sipe al i : dae , 
® 5” % 4 \ 
} ry “| ) Dy i 
u 1 : 2 
Bae ; 4 a ree , oe SOS 
‘i A { sydd th 
ft » : ‘ ry eng ‘ Big 
| ; 
RD eB OLN ere nt raned Wk 
FS et we iy Pa ni i Pee 
Nig ee i RR aie A at Hit 
i if Heri ie 9 i w ) a 
ey tee Ni ud Ae gt SWAG hi Lad ue 
Yiu fi 4 / et Haar ¢ : 4 
‘ ] yh? 4h f a wi Ay Cai Opies abi Ke 4% Hi 
‘ nate Fi 42 ty ; Pats sks Co ts a nae 
As ~ mh Cie ¢ “14% rf noe; wha eis bh ietpink 
ary eee a j Re Sa ai ieee at AGN as 
eae ye ts SERRE DU Int 
Le ri at Od ght Mia ve rain 











) es ae Oo hey) ITE 
7 i ; am, ; mn a i] wig n a 
. 7 a ‘ 


fl Tis 
V #6 














4 
% CONVERGING! 





PrReasoURs ON THE LNTERFACE BETWEEN A 
PVtbeewe DGe NN Deky PAST SE LUT D BOTTOM 
: 
Joe 
| Masami Kawamura ) 
and 


roanaALTsS LOannow 


December 1978 





Mer? SOT, & y BR 


Poogev-amsetpuolac release- distr2bucion unlimited. 


11874 


may if \ 
, 


ren 
J wy, 





ENCL ASS IE IED 


SECURITY CLASSIFICATION OF THIS PAGE (When Dara Entered) 


REPORT DOCUMENTATION PAGE 


. REPORT NUMBER 2. GOWT ACCESSION NO. 


4 TITLE (and Subtitie) 













READ INSTRUCTIONS | 
BEFORE COMPLETING FORK 


RECIPIENT'S CATALOG NUMBER 





a 








3S. TYPE OF REPORT & PERIOO COVEREO 
Master's Thesis; 
December 1978 


. PERFORMING ORG. REPORT NUMBER 












Pressure on the Interface Between a 
Converging Fluid Wedge and a Fast Fluid 
Bo COM 











8. CONTRACT OR GRANT NUMSERa) 





7. 





AUTHOR( 8) 


Masaml Kawamura 
jToannis hoawine U 







10. PROGRAM ELEMENT, PROJECT, TASK 


ATION NAME AND ADORMESS 
od AREA & WORK UNIT NUMBZAS 





9. PERFORMING ORGANIZ 






Naval Postgraduate School 
Memtkerer, California 93940 






11. CONTROLLING OF FICE NAME ANO ADORESS 12. REPORT CATE = 
A x December 1978 
Naval Postgraduate School 


ea ury © t= : 13. jmaQ BEROF PAGES 
Monterey, California 93940 MOD . ; 
4. MONITORING AGENCY NAME &@ ACORESS(!/ differem? from Controlling Office) 1S. SECURITY CLASS, (ofl thle report) 





lSa. DECLASSIFICATION/ OOWNGRADING : 
SCHEDULE . 


16. CISTRIBUTION STATEMENT (ofl thia Report) 


Pee mevee £LOr public release; distribution unlimited. 


17. CISTRIBUTION STATEMENT (of the ebatract mitared in Black 20, if different fram Report) 


18. SUPPLEMENTARY NOTES 





119. KEY WORDS (Con:sinue on severee side ij necessary and identity ay bia0cz numéer) 


ressure Amplitude 
Pnase Distribution 


20. ABSTRACT (Continue an severse side if neceseary and identify dy adlock muamser) 


The pressure amplitude and phase distribution along the 
interface between a tapered fluid layer and an underlying fast 
fluid bottom were investigated both theoretically and 
experimentally. 

Two different tneoretical models were compared experi- 
mentally: a simple model based on a combination of normal 
Modes and ray theory and an exact solution based on the method 


DD , are 1473 EDITION CF I NOV 6318 OBSOLETE UNCLASSIFIED 


(Page 1) “i : I! SFCUMITY CLASSIFICATION OF THIS PAGE (When Data Irivsrec) 





UNCLASSIFIED 


SECUMITY CLASSIFICATION OF THIS PAGESWhen Nore Entared 





mee ABSTRACT Cont'd) 


of images. The experiment was conducted at 100 kHz with a 
wedge of siliccn oil separated from a large tank of fresh 
water by a thin Mylar diaphragm. The simple model failed 

to predict adequately the pressure amplitude and phase along 
the interface. The method of images gave accurate predictions. 


OD Form 1473 = 
Mian. 7 eee Oe 


ee ae i gh i: ae 





Approved for public release; distribution unlimited. 


Pressure on the Interface Between a Converging 
Fluid Wedge and a Fast Fluid Bottom 


by 
Masami Kawamura 
Lieutenant, Japan Maritime Self Defense Force 
B.S., Naval Postgraduate School, 1978 
and 
Ioannis Ioannou 


Lieutenant Commander, Hellenic Navy 
peoey etc cere scoractate School, 1978 


Sepmmeeea im Dartial fulfillment of the 
requirements for the degree of 


Vota oe scl PNCE EN ENGINEERING ACOUSTICS 


from the 


Nave, POSTGRADUATE SCHOOL 
December 1978 





ABSTRACT 


The pressure amplitude and phase distribution along 
the interface between a tapered fluid layer and an under- 
lying fast fluid bottom were investigated both theoretically 
and experimentally. 

Two different theoretical models were compared experi- 
mentally: a simple model based on a combination of normal 
modes and ray theory and an exact solution based on the 
method of images. The experiment was conducted at 100 KHz 
with awedge of silicon o11 separated from a large tank of 
fresh water by a thin mylar diaphragm. The simple model 
failed to predict adequately the pressure ampiitude and 
phase along the interface. The method of images gave 


Meeurate predictions. 
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ie) LINER OeDUCTION 


As sound from the open ocean propagates up onto the 
continental shelf, lower frequencies become lost to the 
water column; energy iS transmitted into the bottom. If 
the bottom has a speed of sound in excess of that of the 
water, then it is predicted [Ref. 1,2] that a highly 
collimated beam will form in the bottom. 

For a channel of constant depth overlying fast bottom, 
normal-mode theory shows that if the frequency of the sound 
1s less than the cutoff frequency of the mode energy can 
be transferred into the bottom. This can be seen if the 
normal mode is represented by upward-going and downward- 
going rays in the layer. As the frequency of the signal is 
decreased, tne incidence (grazing) angle between the ray 
and the bottom increase. Energy continues to be perfectly 
reflected and remains trapped in the layer with no transfer 
into the bottom, until the incidence angle becomes equal to 
the critical angle. At that frequency the mode stops propa- 
gating in the layer and energy is transferred into the 
bottom. This concept can qualitatively be applied to a 
tapered layer over a fast bottom. As a propagating mode 
travels up the wedde, the incidence angle of the equivalent 
ray increases with each bounce until it reaches the critical 
angle, at which point energy will be transmitted into the 


meetom. 





Theoretical and experimental studies of this phenomenon 
have been carried out at the NPS for three years encom- 
passing three generations of experimental setups. 

Lt. Edwards [Ref. 1] with his thesis adviser, Professor 
mee. Coppens Carried out initial investigations in a large 
meee cCliled with brine solution over which a thin layer of 
Olive oil was carefully added. Since the interface between 
the two fluids was horizontal, a wedge-shaped upper layer 
was created by inserting a board of pressure-release material 
at an angle to the horizontal. The sound field was observed 
in the bottom fluid at various distances and depths. When 
plotted on properly chosen coordinates, the data were 
observed to cluster about a Single curve consistent with 
the existence of narrow beam of sound enetering the bottom 
at a relatively shallow inclination. 

In the meantime, Professor Coppens developed a simple 
theoretical model that predicted the shape of the beam 
meccern. 

In an attempt to obtain more precise data, Lt. Netzorg 
met. 2), with his thesis adviser, Professor J.V. Sanders, 
redesigned the apparatus to consist of a wedge of water 
separated from a small tank of brine by a l-mil Myiar 
diaphragm stretched tightly within a frame. The bottom anda 
Sides of the tank were lined with absorbing material to 
reduce reflections. The results showed that the entry point 


of the beam was close to where it was predicted by the 





Simple model. Agreement between experiment and theory was 
fair but not good enough to verify the model. 

Subsequent to these first experiments a more refined 
model has been constructed, with silicon oil in the wedge 
and a large water tank for the bottom. But the results 
obtained in this apparatus by Professor Sanders were still 
Significantly disturbed by what appeared to be some sort 
of interference. 

Experimental results from the three generations are 
compared to the predictions of the simple model in Table I-1l 
Me@erig. I~l. Significant differences between theory and 
experiment were observed on the lower side of the beam 
Dattern. In one extreme example, two veaks of beam were 
observed in a third generation of experiments. Generally, 
the observed angle of depression of the beam is shallower 
and the beam width narrower than predicted by the simple 
model. 

Memciis POlint, the research effort was redirected to 
go back to examining the amplitude and phase distribution 
along the bottom; these quantities provide the link between 
mea] acoustic field in the layer and the formation of the 
beam in the bottom. 

The objective of this research reported herein are to: 

1) Develop a simple, closed form expression for the 
pressure distribution along the bottom that can be 
used to reliably predict the beam pattern to the 


ico ete) ml 


0 





2) Calculate the pressure (including phase) along the 
bottom predicted by application of the method of 
images and compare with experimental results. 

In this report, two different coordinate systems are uSed. 
The first is a special coordinate system applicable to the 
Simple model; the origin is located at a point on the bottom 
where the lowest mode reaches cutoff. The second is a 
general coordinate system more Suitable for the image 
calculations and experiment; the origin is located at the 


apex. These two coordinate systems are shown in Fig. I-2a,b. 
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a ee nO, 


fee OU RE DISTRIBUTION PREDICTED BY A SIMPLE MODEL 
meeeeesic Concept 

The purpose of this part is to summarize the 
theoretical background necessary to conceptualize the 
propagation of sound in a tapered fluid layer above a fast 
bottom, using a combination of normal modes and ray theory. 

a. Layer of Uniform Depth 

One-dimensional propagation in a fluid layer 


of uniform depth H with speed of sound c, and density Oy 


iL 
overlying a second layer with a faster speed of sound 
C5, and density P95 can be exactly represented by an equation 


@e the form 


PB =) Ae Sele 5 woe = (II-1) 
~ n 
n 
Mtece Z' is the depth, x' is the horizontal distance, 
Ko and Kos are the components of the propagation vector 


Ll 1) 


— 
fror the ace mode in the z' and x' direction respectively. 


AL is the magnitude of the not mode. For propagating modes, 

the propagation vector K can be represented as in Fig. II-l. 
es 

The angle between k and the horizontal is @. Snell's law 


gives the relationship between the speeds of sound in the 


Mayer and the critical angle oe 


IS 








— 


merge Li-l. The propagation vector k and its components 


y) = arc cos = (II-2) 


As long as 6 < oF the energy is totally reflected. How- 
fo, When @ > 9. energy is transmitted through the 
previously totally-reflecting surface. 

To satisfy the requirements of propagation, 


Ko must be a constant, real number for each mode and from 





Ber. 3, 
Ko 
a ea 
fam Kot H = ~Alrmaed | ) (iLi—3) 
n a3 YH 
where 
2 ls 7 
Cee x! ‘So? a” 
2 
me Cutofr @ = DG and then the geometry of Fig. II-l 
implies k_, = ~~ so that 
Me C5 
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Cope nee — hy = =) 7 (II-5) 


and for a mode with frequency just above its cutoff 
frequency EW there 1s a pressure antinode at the bottom. 


Now by writing Eq. II-l in the form 


A Eh era ~1ik,, rhe! 
a ) a be n 25 n y qi (wt-kx ) 
n 
oc Geaex tk, 2!) 
a = ) iA fe “a ee 
2 nt 
n 
ie(itemica a sc Ks 
eet are 
-e } (6s) 


we have two plane waves that propagate with wave number 
~ 
k = kx +k ,2'. If we assume that the frequency is 
n rel 
Buewly lowered, higher modes are cutoff as the frequency 


decreases with Kot sO SiglereKer ebL ele Oe yion =) 
n 


b. A Tapered Layer 
Let uS now consider a wedge with angle 8 as 
indicated in Fig. II-2 and II~-3, where x 1s the distance 
from 9 = ore to the apex, and xe is the distance of closest 
meeroacn (Turnaround distance) of the indicated ray. In 
mms Case each bounce off the bottom increases 9 by 28. At 
the point where $ = OO energy begins to penetrate into the 


lower layer. The resulting attenuated ray bounces off the 


ey, 





IP oblope 


Ti~-2. 


. 


The geometry of a wedge 
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top and makes another collision with the bottom, losing 
still more energy. This process continues until the wave 


reaches the "Turnaround distance" Xe defined by: 


x, = Xo(i- sin 9) (Ty 


For the wedge then, the frequency remains constant while 


86 increases with each bounce off the bottom until that depth 


A 1s reached for which 6 = oF where 
(n-=x=)T 
Ay = Se (TI-8) 
id 
n 


meem Fig. I[I-l we obtain at cutoff 


where 


kK = & (II-10) 


M@emesuostieuting Eq. (I1I-9) into Eq. (II-8) yields 


(n-5)A 
— WP es 
aD Cane, a 
es 
Knowing Hy allows calculation of the value of x' = Ko where 


energy will begin to be transmitted into the bottom. 
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2. Pressure Distribution along the Bottom 
A.B. Coppens [Ref. 4] formulated an approximate 
model to predict the characteristics of the transmitted 
beam in the bottom. The expressions for the pressure ampli- 


tude and Kye are: 


A = Ag S Cie eZ } 
ren = = = ee = 
Viole 





where x', Xa are defined from the geometry of Fig. II-l, 
9 meom Eq. Li-2, A, 1s the pressure amplitude at x' = 0, 


and A defined as 











7 E 1 OVA ae 
i = ; a fare sin OSes arc sin Ter 
b°-1 
IL Je Cat OZ One 
io Sar ae sin o77z) Bligeh SeGhiei men) 
l-c 
i. = e eas 
+ 3 Bei stare Ha agegy > Oe Genres 
le (TI-14) 
where 
a (21-15) 
Py 
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Q 
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Q 
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CD 
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(ii. 


Z = l1- (1 - =) TE esl) 


In the case where 


meecq further 


pee Zemin a 





» 1 Ss Be es 
P45 8 3 Xo 
A = A, e (II-18) 
and 
C =e Kae ks) = =x 
x 
Sane 
9 = arc sin (=) (e209) 
Le Sy 
x 
O 
3. Pressure Distribution under the Bottom 


In this part we present without detailed derivation 
A.B. Coppens [Ref. 2] work concerning the pressure distri- 
bution under the bottom based on the simple model and 
obtained through the use of the Green's function and the 
method of stationary phase. Two basic assumptions were 
considered for this development. 

a. The plane wave has wavefronts that are parallel 


to the shore line. 


PPA 





b. The pressure distribution on the bottom can be 
Viewed as an acoustical source which generates 
the beam of sound projected into the bottom. 
The geometry 1S indicated in Fig. II-4. 
The Green's function formulation for this case can be 


developed through the geometry of Fig. II-4 as follows: 


At the point (x',z') the instantaneous pressure is given by 
P= Bea aw Ge) Ge as (II-21) 
where 
ies ny ee dx" 2 
pest) = ef Byte Srl sag (11-22) 
The integrand of Eq. (II-22) is formed from the pressure 


meestripution at the bottom 





Ceara (it\=24) 


Zo 








PLANE OF SEA SURFACE 
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Fig. II-4. Plane wave propagation in the wedge 
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imeeeg. (LI-23) the term _—s ioc cn moO come Correct ion 


x 
Ll - yr 
O 

for the cylindrical convergence of sound. In the absence 


of losses into the bottom and if there is no energy loss 
within for the wedge, the pressure amplitude would behave 


as 


ie 
2+ = = (II-25) 
H 
O 
also from the geometry of Fig. II-2 we have 
H 
a sr (II-26) 
oS ran 
O 
q 
so that 
e 
H | om 
2) i - xX 
yi 
The term A(x") is the pressure amplitude on the bottom 


and the exponential term is the corresponding phase. 
This pressure distribution is an integral of the 


general form 


dx (in 72 oy) 
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where F(X) varies slowly compared with the exponential 
Factor. The value of the integral can be well approached 
with the method of stationary phase. 

Using the above method and assuming the simple 


u 


model for oT << = cos*é., and for small 8 the following 
O 


expression for the pressure distribution was developed: 





a L 3 
oe 
\P| = oe or we - 2c (II-29) 
. Py O 
where 
a WCU cos Q ) Sin g (II-30) 
ean 0 teal oie 
eee eressure Distribution along the Bottom Predicted 


by the Method of Images 

Since the simple model did not yield results agreeing 
with the experimental measurements, it was decided to use 
the exact but more cumbersome method of images to attempt 
Mempredict the pressure distribution on the bottom. 

Let the source be a distance X from the apex of 
the wedge and a depth H below the sea surface as indicated 
meraig., Ll-5. 

Define 


arc sin CT 37} 


ic 


mo 


4 
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and 


o, 2 B- 4 (I-32) 
where 8 is the angle between the sea surface and bottom, 
bg is the angie formed at the apex between the source and 
the sea surface, and 6 1s the analogous angle for the source 
and the bottom. 

Hits on is the angle formed at the apex between the 


nth image and the bottom then Fig. II-5 reveals 


oe - ae oo = 2B = $, 
eo yk 2(8-9.) ee oe 
> = 36, + 4(B-o)) = eave (i 2 
Sesmore generally 
o, = (ntlo, + n(8-o5) = mB + o> for n even 
oe = ere (nt+1) (8-9) = (nt+1)6 - We omen .odd 


l= 34) 
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Fig. II-5. Geometry of a wedge by the method of images 
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This set of equations can be reduced further to one 





expression: 
m. Ie ae haen = elas IL Ea 
@, = int 5ll-(-1) 338+ (-1)"o, = 2 INTI==18+(-1)"9, 
(II-35a) 
Also for the symuetric member n' of the paired images 
we obtain 
= texas ay ge gl = oa ile ae, - Dee ons ole 
at {n eye (a ae Care) eS Ae tUlae | 5 oe ee) os 
ILE Ss Bey, 
where INT[ ] denotes the largest integer which is equal 


Ee, Or Smaller than the argument. 

Considering the geometry of Fig. II-5 and II-6 we 
can derive expressions for the distances rn and ar which 
represents the distance between a receiving point on the 

ele ee! 
bottom at distance x from the apex and the n and n 
image respectively. Also we can derive expressions for 
the angles a and is which formed at the receiving point 


between the bottom and the images n and n' respectively. 


These expressions are: 


r, = xX" + i - 2Xx cOs®., (Palle: 


Zo 











Fig. II-6. Geometric development of the image's position 
with respect to the receiving pein e 
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Oy, = eer can a's (iat= 37-2) 
n x 
Sem 5) 
6, = arc tan [ j (Cio) 
ae cos 9 ae 
Ey x 
Let us now define 9 euicm@on, . as the angles of 
nm me an 
incidence for the qe bounce from the bottom for the n 
and n' image respectively, m= 1,2,3, ... (The ate bounce 
is the last-one before reaching the receiver.) From the 


geometry of Fig. II-5 and Fig. II-7 we obtain expressions 


mer oO as follows: 
nm 


Ta, 7 ay 
7) een) aa 
yee on 26 Grea a) 
102 Sa 
ceMMER 5G oo” 


The general expression is 


; Ss eer (II-39a) 
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II-7. Geometric development of reflection angles 


ge 





Also for the symmetric member n' of the ae pair we obtain 


9 = 6 Se Zauils i= 3.91) 


and the maximum number of bounces of the ae and ee 
image is 

4 o8 ee ale eee 

Bee NE a a ae 
a ue 1 a - 16 es ad 
= 5 (n tL as Z L -{ (-L) } 
oe oes 
= INT i53! = INT l=] (in 0) 


[mene receiving point is on the bottom, then from Fig. II-8 


we obtain: 


1. The distance from each member of a pair of 
images (n,n') to the receiving point is the 
Same selec ec. =... pe 

n n 

2. The number of reflections (from the bottom) 

from the Ae image to the receiving point 1s 


One more than the ae image. 

3. The grazing angle to the bottom is the same 
for each corresponding bottom bounce so that 
the reflection coefficients for corresponding 
bounces are equal, except the last bounce for 


the nue image 


So 








“ = 
a “ = 
, \ 
“XN 
‘\ 
x 
» 
i \ 
\ 
. 
\ 
\ 
D PN \ 
\ 
\ 
AO 
; | 
! 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
A 
4 
4 
YA 
7 
we 
a“ 
“ 
Zz <a we 
=e —eet 0 sel ry’ 
\2 
Fig. II-8. Geometry of symmetric images 
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Consequently R = R_,._ where R and R_,. are the reflec- 
~nm ~n’ m ~nm ~h Mm 











tions coefficients for the ma bounce of the ee and i as 
image respectively. 
The reflection coefficients for the ea image are 
given by 
OxC 
DG y 
PCy ~nm 
“nm pc, (II-41) 
LoD 
Pio] ~nm 
where 
e 
EE (<4) cos*6_ 
ee eee 6«6Cfor 6 Cl 8 (II-42) 
Y = Sin 6 nm — “ec 
~nm nm 
and 
v = Oi < 6 
va, nm C 
(tT=—43) 
with 
el 
) = arc cos (—) (II-44) 
e C5 
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Let uS now consider the pressure P(x) at a point x on the 


bottom. The contributions from the paired images are 


N 
ecicce int (2S 
ice ee) } (II~45) 
mee One =] 
and 
N ] 
oe Int (2) 
) {—e T Ra (-1) } (II-46) 
re ~ Tl 
i Onna n=0 
(The nth and nth images differ only in that Ae image has 


no R reflection.) 
Pt1O 


The complex pressure 1s 


~no m=0- Premera) 


A computer program was constructed to calculate this 
complex pressure (amplitude and phase). This program, 
with the corresponding flow chart, 1s presented in Appendix A. 
imeecase OL aesource SO Gistant that it can be 
assumed to be at r +> ~ then the geometry of Fig. II-9 


yields the approximation 
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Fig. II-9. The geometry of distance source -hojebaop aeit- gen mere! 


aif) 





s) ae) (II-48) 
which simplify the mathematics Significantly. Here 


Rio NOE = R(9 ) (Moths) 


In this case, the expression for reflection coefficients 


are (see Appendix B) 


M 
gn om ~ ~O ~ R$) 
M 
Ri i: ek 7 Rt, ) 
Sh 
TE Bom 7 Rok = Rlog) RC) (II-50) 
M 
eames, pty RAO 5! 
M M M 
ep enm ~ gan 2m an ~ gens 2m R(O,) 


Papere the right hand sides of Eqs. (11-50) depend 


only on n, then define 
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12 Cpa ge (0 
R ime Me Tl Chis i.) 


[(n-2)R, POG ete 2 ac 
Also, the geometry of Fig. II-9 gives 
1g = r-=- X cos ae (TI-52) 


Substitution of Eq. ((II-50), (II-51), (II-52) into 


mae {LI-47) yields 


ikxcos?, INT [———] 
iE 


The last expression for pressure is the approximation for 
a distant source, which is described in detail elsewhere 
beet. 4). 

iewadvantage Of Eq. (11-53) with respect to Eq. (II-47) 
is that because the elimination of product term, the 


calculations are greatly simplified. 
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Initially it was intended that the beam pattern in the 
bottom would be exhaustively investigated, and a larger 
experimental system with silicon oil in the wedge and a 
large water tank for the bottom was deSigned. But, even 
in the new system, 1t was observed that the beam pattern 
was being significantly disturbed by unknown interferences 
that caused the apparent angle of the beam to move to 
Shallower angles than predicted by the simple theory des- 
cribed in II-A and made the beam seem narrower than predicted. 
This interference was definitely not from the walls or 
bottom of the water tank. It might have been associated 
with the wedge, but whether it was an inherent property of 
wedge propagation omitted from the theory II-A or merely 
an artifact of the experimental setup was the question. 

In an attempt to study this question more thoroughly, 
the experiment was redirected to investigate the amplitude 


and phase of pressure on the bottom of the wedge. 


A. APPARATUS 

The experimental facility consisted of a wedge of sili- 
con 011 overlying a large tank of fresh water; the two 
Fluids being kept apart by a thin Mylar membrane. A photo- 
graph of the apparatus is shown in Fig. III-l. The sound 


source, positioned at the end of the silicon-oil wedge, 


sent a pulsed Signal toward the apex of the wedge. It was 
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possible to position the receiver anywhere within the fresh 
water region, but in this experiment the receiver was mainly 
placed close to the bottom of the wedge. The geometry is 
Swewoaein Fig. III~-2. 
me Water Tank (Fig. III-3a,b,c) 

Tne water tank, constructed of 3/4 in. thick 
plywood, had the inside surfaces painted with Varathane 
but were otherwise uncoated. It was filled with tap water 
at room temperature. 

2. Wedge (Fig. III-4) 

The silicon oil was contained in a frame whose 
Sides were made of 1/2 in. x 2 in. aluminum and whose 
bottom consisted of a thin Mylar sheet. Half-mil thick, 
aluminum-coated Mylar was fastened to the bottom of the 
frame by first taping a sheet of Mylar, metalized side up 
to a flat surface. Then the aluminum frame was placed on 
the Mylar, and the edges outlined by a Sign pen. A bead 
of "Royal bond GRIP" contact cement was applied to the 
bottom of the frame and to the marked areas of Mylar. frer 
the cement was completely dry (about 15 minutes), the 
frame was placed on the Mylar and about 50 pounds of weight 
Beaced On top of the frame. After the contact cement had 
Mec eor '2 nours, a coat of “Silver print conductive paint" 
waS applied to the outSide joint between the Mylar and the 
frame. This provides electrical contact between the Mylar 
and the frame. (This was necessary to allow measurement of 


the slope of the wedge as will be explained ina later 
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section.) The excess Mylar was trimmed flush with the 
outside edge of the frame by a razor blade. As a final 
step, a thin coat of the "Silicone rubber sealant" was 
applied to both the outside and inside of the Mylar-frame 
interface to prevent leaks. Because the tension in the 
Mylar tended to bend the walls of the frame inward, a 
3/8 in. diameter aluminum rod was used as a "frame spreader" 
to spread the frame and keep Mylar flat. The length of 
the frame spreader was adjustable. It was placed across 
the shorter dimension of the frame and adjusted until the 
Mylar film was mirror flat. Aluminum arms at each end were 
used to suspend the wedge from the top of the water tank. 
The longer arm was connected to the apex side of the wedge 
and extended in the propagated direction. The slope of the 
wedge was set by raiSing the end of this arm with wooden 
blocks. The other arm, perpendicular to the propagation 
direction, was connected to the source side of the wedge. 
The horizontal level of the wedge was adjusted by means of 
screws at the ends of this arm. 
be cOUrce Transducer 

The source transducer was made by Lt. Netzorg. 
Construction details are discussed in his thesis [Ref. 2]. 
The design of the transducer was quite well done: it has 
Peericient horizontal directivity so that sound is not 
reflected from the side of the frame, there is negligibie 


sound propagation in the backward direction, and its active 
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face, 2 cm high by 8 cm wide, is the largest size allowable 
in the wedge. The experimentally measured directivity is 
very close to the directivity predicted for a rectangular 
Beumeece (Fig. [ili-5a,b,c). The theoretical directivity 


is described as follows: 


sin (5kL,siné) sin (SkL,sind) 

fio ,d) = =. = (Ref. 3,5] 

ZKL, sing xKL,Sind 

(tees = )])) 
where 

L, = Horizontal dimension of an active face 

of a rectangular source. 
L, = Vertical dimension of an active face of 


a rectangular source. 


Ww) aT 
K = = = ae Wave number 


Q = Horizontal angle measured from a major 
radiation axis. 


0) = Vertical angle measured from a major 
radiation axis. 

As a matter of convenience, the directivity measurements 
were done in the water instead of silicon oil. The fre- 
gGuency used in these meaSurements was adjusted to give the 
Same wavelength as produced when the source operates at its 
prescribed frequency in silicon oil. A frequency of 
127.8 kHz in fresh water (Cc. = 1503.8 m/sec) produces the 


same wavelength as 100 xHz in a silicon oil (cy = 1176.5 m/sec). 
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For more accurate comparison between the experimental 
results and the method of images, the vertical directivity 
was employed in the computer model "WEDGE O" (APPENDIX-A). 

The source transducer was mounted on a support that 
allowed it to be located anywhere in the upper fluid layer. 
Vertical adjustments were made by means of a micrometer with 
Mmeeein. Of travel. 

4. Receiver and Attachment (Pig. III-6a,b) 

A LC-10 hydrophone, manufactured by Celesco Indus- 
tries, was used as the receiver. It was possible to locate 
the receiver anywhere in the fresh water. Since the main 
purpose of the experiment was to investigate the sound 
pressure distribution on the bottom, the hydrophone was 
usually located as close to the bottom as possible. The 
attachment mounted on the frame of the wedge as shown in 
mog., ITli-6a. 

At the beginning of this thesis, some measurements 
were made of the beam pattern in the lower fluid. Then a 
different attachment was used (Fig. III-6b) which was mounted 
on the large arm. 

In both cases, the depth of they hydrophone was 
adjustable. In all cases, measurements were made in the 


Wermcical plane that contained the acoustic axis of the source 


transducer. 
5. Slope-measurement Assembly (Fig. III-7a,b,c) 


The slope measurement assembly consisted of a pin, 


micrometer, and attachment. Fine vertical adjustments of 
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the pin were made by means of the micrometer with 1/2-in. 
of travel. The attachment was mounted on the wedge frame 
and could be located anywhere in the wedge. The pin was 
electrically insulated from the frame, so that when an 
Ohmmeter was attached between the frame and the pin, the 
precise depth at which the pin contacted the metalized 
Mylar could be determined. 
See Electrical Setup 

Figure III-8 is a block diagram of the electrical 
setup which consists of a driving system, a receiving system, 
and a measuring system. The transducer was driven by a 
General Radio Oscillator (Type 1310) and Tone Burst Generator 
miyee 1396-A), polarized by a 300 Volts D.C. biasing-voltage 
[Ref. 2]. The driving frequency (100 kHz) was measured with 
a Hewlett Packard electronic frequency counter (Type 5233L) 
and the driving waveform displayed on one channel of a 
Tektronic Dual-Beam Oscilloscope (Type 565). The received 
Signal was amplified (40 dB) by a Hewlett Packard Amplifier 
(Type 465-A), and band passed (40kHz ~ 200 kHz) through a 
Spencer-Kennedy Laboratories Band Pass Filter (Serial 1683 
Model 302). The received signal was displayed on the other 
channel of the Tektronix Dual-Beam Oscilloscope (Type 565) 
and the amplitude was measured on this display. The received 
Signal also went to a Hewlett Packard Oscilloscope (Model 
120B) where its phase was compared to that of the driving 


Signal. 
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fe AUxtliary Apparatus 
To measure the sveeds of sound and densities of 
fluids, other pieces of equipment were used. 
Transducer 
Two Hydrophones LC5-2 
Analytical Balance 


ASTM Vessel 


Pee — ROCEDURE 
Before beginning the experiments, the speeds of sound 
pm@emwcaensities of the fluids (a silicon oil for the upper 
layer, fresh water for the bottom) and the slope of the 
bottom was measured with high accuracy. 
1. Sound Speed Measurements (Fig. III-9) 
The speeds of sound were measured in a pan which 
Maa CGCimensions 10 3/4 in. x (mone es oO ein. About 1500 Ze 
of a test fluid was carefully poured into a pan to prevent 
the entrapment of air bubbles. Two LC5-2 receivers and a 
transducer were used to measure the speeds of sound. The 
distance between two receivers was set at exactly 20.0 cm. 
And the transducer was driven by a 150 kHz pulsed signal. 
The time of transit for one receiver to the other was measured 
using the time-delay feature of the dual-beam oscilloscope. 
The speed of sound was then calculated by dividing the 
distance by the transit time. 
2. Density Measurments 
The densities of the two fluids were found by first 


weighing a 100 ml, ASTM vessel on an analytical balance 
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readable to one-tenth mg. The vessel was then filled to 
the 100 ml mark and weighed again. By subtracting the first 
weight from the second weight, the weight of 100 ml of the 
test fluid was determined. The weight of the test fluid in 
grams divided by the 100 ml volume and, the quotient times 
1000, gave the density of the test fluid in kg/m?. 
3. Measurement of the Slope Angle 

The approximate slope was set by raising the end 
of the long arm a precalculated height with the use of 
meeaen blocks. Then 3500 cc of silicon oil were slowly 
poured into the wedge, and the frame-spreader was adjusted 
to keep the bottom flat. The water level in the tank was 
then adjusted to keep the apex line (shore line) straight. 
The source transducer was set at about middle depth on the 
wedge centerline. The slope angle was accurately measured 
with the slope-measurement assembly described in the appara- 
tus section. Because the top and bottom of the frame of 
the wedge are parallel, if the Mylar sheet is flat the 
distance from a reference level to the bottom should be the 
Same at all points. The distance from a reference level to 
the surface and to the bottom were measured at several points 
on the wedge centerline, and the data were plotted. The 
Slope was calculated from the plot as shown in the example 
Mmetustration in Fig. III-10. 

4. Measurement of Pressure Amplitude along the Bottom 
The electrical setup of the measurements was shown 


mieeig. tit—-8. The oscillator output, set at 100 kHz, was 
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fed into the tone burst generator. A pulsed signal of 
36 cycles of the 100 kHz input with a closed-gate duration 
time of 10 msec (repetition rate was about 100 Hz) went 
through the transducer biasing network and then to the 
source transducer. The driving signal was displayed on one 
channel of the dual-beam oscilloscope and adjusted for a 
peak-to-peak voltage of 10 Volts. The receiving signal was 
displayed on the other channel of the same dual-beam oscillo- 
scope. This was very convenient because the receiving sig- 
nal was always compared with the driving signal on the same 
oscilloscope. Reflections of the receiving signal were not 
observed. 
5. Measurement of Pressure Phase along the Bottom 

The phase of the receiving pressure was compared to 
mee of the driving signal on an oscilloscope. Fig. III-8 
also shows the electrical setup for these measurements. The 
Signal from the oscillator was directly fed into the x-axis 
of the oscilloscope. The receiving signal was amplified 
and filtered, then fed into the y-axis of the same oscillo- 
scope. The signal on the x-axis was continuous with a con- 
Stant amplitude (adjusted to 10 Volts peak-to-peak voltage). 
On the other hand, the signal applied to the y-axis was 
pulsed and had an amplitude that depended on the receiver 
position. The pattern observed on the oscilloscope was an 
ellipse crossed by a line on the x-axis (like "$5" of the 
Greek alphabet). The line occurred during the time the 


received pressure was zero, i.e., between pulses. When the 
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Signal, generated by the oscillator, was sent directly to 

the transducer, bypassing the tone burst generator, the 
pattern on the oscilloscope was the same ellipse but without 
the line on the x-axis. This indicates that the pulse is 
sufficiently long to allow accurate phase determinations. 

(It also showed that reflections in the tank were negligible.) 
But measurements were still difficult, because the amplitude 
of the receiving signal was not constant and the phase 
changed very quickly with changes in the receiver position. 
The phase measurements were done independent of the ampli- 
tude measurements. The position of the receiver was recorded 
when the phase was 0°, 90°, 180° and 270°. (Corresponding 

to a straight line with positive slope, a symmetric ellipse 


and a straight line of negative slope.) 
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Fig./l[-5b Horizontal directivity pattern H(@,0) 


and boundary of the wedge 
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Fig.[l[-5c Vertical directivity pattern H(0,9) 
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Fig.[[-6a Receiver and attachment (1) 
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Fig.[ll-6b Receiver and attachment (2) 
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Veen ce eo aN “COMPARTSON 


Measurements were conducted under two different 


Senaitions: 


Common Parameters 


Speed of sound in the wedge oS 
(ilicon) 7 
Speed of sound under the bottom Cc = 
(fresh water) 
Bemsity Of a Silicon oil P, = 
Density of a fresh water ig = 
Frequency f = 
Puese length Lo = 
Closed gate duration te = 
Pulse repetition rate - = 
eAcst L 
BuGpe angle 8 = 
Source distance from the apex X = 
Source depth H = 
CASE 2 
Slope angle g = 
Source distance from the apex Xx = 
Source depth H = 


The speeds of sound and the densities of the 


127625) m/sec 
1503.8 m/sec 


Ooo “g/m 
99 5.::0 kg/m 
100 KHZ 

36 wavelengths 

10 msec 


POO HZ 


75.5 Cm 


io. Cm 


fluids were 


measured at 18.5°C. The two cases differed mainly in the 
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Slope angle: for the first case 8 = 2.7° and for the 
second § = 2.4°. These are fairly realistic values: 
@eeording to Sverdrup-Johnson’s - Fleming [Ref. 6], con- 
tinental slopes off mountainous coasts have, on the 
average, a Slope of about 3.5°, whereas off coasts with 
wide, well~drained coastal plains, the slope is about 2°. 
However, there are many different estimates based on differ- 
ent measurements at different locations. For extreme 
examples, estimates of the average continental slopes range 
[mmem about 5° to about 0.1° [Ref. 7]. 

The distance from the apex to the source was limited 
to less than about 6 times the lowest possible cutoff 


distance Xo (XO =o. 92 em £Or CASE 1, Bee ios tt Or eo i 2). 


esse orn lLMENT VS. SIMPLE MODEL 

1. Pressure 

Figures IV-la,b show the experimentally-obtained 

pressure distribution along the bottom. In both cases 
there was a dominant pressure peak near the lowest cutoff 
point but shifted slightly in the direction away from the 
apex. To compare the simple model with the experiment, 
data and the theoretical curve were plotted on the same 
Peeaphn fOr the interval 0 < x < 2x, (hig.  iv-2ey,5). There 
are Significant differences between the experimental results 
and the simple model. The theoretical curves decay quickly, 
and the model predicts that there should be no energy 


transmitted into the bottom, outside the region are Sig x, 
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(the small amount observed in the simple model results 
from the approximations). On the other hand, the experi- 
mental curves have a long tail with energy transmitted 
continuously into the bottom all the way to the apex, and 
a pressure peak shifted significantly, so that energy is 
leaked into the bottom before the cutoff depth, x > Xo 
2. Phase 

The phase distribution was measured in the region 
O< x < 6x, for CASE 2. The rate of change in the measured 
phase was almost constant but with small irregular fluctua- 
tions. Figure IV-3 compares the experimental results within 
egeeregion 0 < x < 2X, With the results of the simple model 


meeemlated from Eq. (II-13) 





For comparison with the experiment, the phase was restricted 


to the interval -7T to +7. 


Figure IV-4 shows the apparent angle of incidence 





Sin 8 Sob gies, 
§ = arc sin ( =) = arc sin ee 
1 - yr xy 
O ie IDV S16) 
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which is calculated from the phase (—k, x!) as discussed 

in Section II-A2. The graph shows that comparison between 
the simple model and the experiment is better accomplished 
by studying 8. The rate of change in the phase is almost 
constant when x > X)J1 Di astemunerdent angle 9 changes 
Significantly in this region. The incident angle @ increases 
with decreasing distance x. At x = x the incident angle 

8 becomes the critical angle aie and Eq. (IV-1) shows that 
-k_,x' = 0. Subsequently the incident angle increases more 
rapidly as x decreases from a ie) Xo° Seecnie Lurniasound 
point x = a the angle 8 becomes 90° and again “kX! = 0. 
Thus, energy can be transmitted into the bottom only for 
oe Xx < Xx, Hee Ggerictmee Lie Model =m PEeecltse Experimental 
investigation of the incident angle was difficult because 
the rate of change in the phase was very small. Significant 
differences between experiment and the simple model were 
observed. From the experiment, there is no turnaround 
point. The experimental pressure distribution has a long 
tall which means energy can transmit into the bottom all 
the way to the apex. And the slopes of the phase curves 
(Fig. IV-3) show significant differences between experiment 


and the simple model. 


pee SAPERIMENT VS. METHOD OF IMAGES 

The computer program "WEDGE O" was designed for a com- 
parison between experiment and the method of images. The 
program has the flexibility of an adjustable source distance 


and inclusion of the vertical directivity of the source. 
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1. Pressure 

Before beginning the comparisons, the computed 
pressure distributions were observed with several combina- 
eeems OL distance (x = 100 m) or near source (x = 1m) and 
with directional or omni-directional source. In all cases 
the major pressure peaks and tails were similar; the 
distance-source case showed slightly smoother pressure 
distributions and was little effected by source directivity. 
On the other hand, directivity becomes important if the 
source was near. As a result, computations were done with 
fPeemaccual Source distance (x = 70.0 cm for CASE i, 
jo 75-.5 Cm for CASE 2) and source directivity described 
in Section III-A3, and these results were compared with 
the experiment. 

Figures IV-5a,b show the comparison between experi- 
ment and calculation. Both pressure amplitudes were normalized 
at the maximum pressure peak. In both cases, there are some 
quantitative differences, but good qualitative agreement; 
most of the local maxima and minima occur at the vroper 
positions and with the approximately correct amplitudes. 

Figures IV-6a,b are enlarged graphs of Fig. IV-5a,b 
memethe region 0 < x < 2X0: From these graphs, it is 
obvious that the tails of major peaks of the experiment 


and theory agreed very well. 
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2. Phase 
Theoretical calculations were also done with the 

computer program "WEDGEO" for a distant source. The 
results were compared with experimental data (Fig. IV-7, 
CASE 2 only) with both phases normalized at x = ee an 
the region oe SxS 2X01 these phases were almost identical, 
but some differences were observed for x < Xi the experi- 
mentally observed phase changed somewnat more rapidly 
than that of the theoretical prediction, the difference 
increasing as the apex was approached. There are some 
possible explanations for this. The depth of the wedge 
at x = Ko was 4.7 mm, and most differences were observed 
Memes the depth was 3.0 mm or less (0 < x < 6 cm). In 
this region, the pressure amplitude was very small 
(Fig. IV-2a), so that measurement of the phase was quite 
difficult. Also, very close to the apex the surface of 
the silicon oil might be affected by surface tension and 
not remain ideally flat. 

The incident angles $ in the upper layer were approxi- 


mated as follows; from Fig. II-l, the angle of incidence 





Ss , 
Koo 
ae a ee eCos Coxe 
where 
k= 2 = 
1 C 





q ! 
Two adjacent points x, , X5 for which KX! have the 
Same values in the interval -1m to +m yield a change in 


phase of 27, 


me cesta ee Xy = 27 
2 2 a 
Ky 7 ors oe oS 
Ll 2 
and 
=a 
0 = arc cos (za) CEVv=2) 


Figures IV-3a,b show the results for the experiment 
compared to the results obtained from the method of 

images and from the simple model. There is no turnaround 
point (for which we would have § = 90°). In both cases 
the incident angle equals the critical angle a little 
further away from the apex than x/KO = 1, which is consis- 
tent with the measured behavior of the amplitude. For 

oo > Xs most of these incident angles were less than the 
@eitical angle re This does not imply that there is no 
energy transmitted into the bottom: these incident angles 
result from the phase-coherent interferences of a large 
number of images, the fields of which have angles of 
incidence on the bottom encompassing a range of values, 
some of which exceed 2a This shows a significant differ- 


ence between the method of images and the simple model for 


Teak 





which energy can not be transmitted into the bottom when 
the incident angle is less than the critical angle.) 

The phase as predicted by the method of images and 
the experimental phase is nearly linear in x. The ampli- 
tude varies more slowly, but perhaps not slowly enough 


for the method of stationary phase to be applicable. 


See aL lEMPTED RECONCILIATION 
The predictions of the simple model do not agree well 
with the experimental measurements. However, investigations 

of the pressure field in the bottom would be facilitated 
if a simple analytical expression for the pressure distri- 
bution along the wedge-bottom interface could be formulated. 

Consequently, the purpose of this section is to com- 
pare the simple model with the method of images for those 
cases where the source 1s very distant. The amplitude OF 
the pressure distribution predicted by the simple model 
mes the form 

OK tand . k as 


Dale 2 os 
a) ( 3 ) "Kk, (yy) 


A = A, e Ci ==) 





where the predicted values are Ky - kK =1, k, = —=— and 


ae , , 
Ky = >: Changing Py7 P51 Cyr Cy and b it was mee 


the pressure amplitude depended only on the ratios = ; 
tan 9 2 
ES 
8 
Next comparing the predictions of the method of images 


with the simple model also showed that the value K4 aaa) 


Cc 


2 


Ve 





gave the proper spatial behavior, to reasonable accuracy. 

This was accomplished by taking ln of both sides of 

Bee (t¥-3) and plotting then ln “0 Vs er 
Obviously the slope of the Bion 


On log-log paper. 





A, Oy Ky tan a Ky e 
a. = as Oe (nl on s } ~ em 
intin() } Ini (5 ) ( Z ) K3} ms Ky In(y ne 
Z Oo 
(IV-4) 
SLvVes Ky: 
; 
By the same technique the exponents k, and k, of — 
tan 6 t : P2 
and —y— and Ene  -Eactcor K were examined but in these 
cases the results for Ky Kos K were not constants; but 
| 2a pial x! 
functions apparently dependent on —, , 3, ore 
Po 2 “oO 


We can therefore conclude that the simple model is 


not generally valid for arbitrary values of Cyr Par Cys 


a 


Ke 
O 


K, ir the simple model predictions are brought into good 


agreement with the experimental measurements for the values 


8 and However, it is useful to note that by assuming 


Se Oils gr/om> 


oF = 

_ 3 
P5 = Ceo Js7gr/ ecm 
e277 ie oem) SEC 
cy = S50 376 7 SEC 


WS 





Dee OO 6 ead 


After the discovery of this agreement an effort was made 


to specify the ranges of variables in which the amplitude 


prediction of the simple model is valid with K 3 = = 
The simple model requires 
oe ih Zz 
44 > ee > Cas i De) 


Assuming that the amplitude distribution is unimportant 


1 
= Such that = ee =, this inequality must 
O O 


be achieved within the restriction (IV-5). 


for those values 


Py tan 9 l 5 
an an a= COS § 
As gn * Bv2 26 Qt (IV-6) 
A 
O 
or 
e) 
Baa < i Sin 6 (1 - ataee ) _* (IV-7) 
BD) . See 


This inequality was tested on the interval 0.002 < 8 < 0.005 
Zor on = 60° and = 0.5 or 1.0 there was uniformly good 
agreement between the simple model and the method of 

images for = =O. ouee TOT = = 1.0 there was strong 
disagreement between the two models. 
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a | 


Therefore the range of validity of the simple model 
with kK = 7 appears badly limited. Further, since the 
phase of the simple model does not agree with the experi- 
mental measurements and the predictions of the method of 
images (as already has been mentioned in part IV-A2), 
we can conclude finally that this model is very weak and 


its application very restricted. 
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V. CONCLUSIONS 


Two main conclusions should be emphasized. 

(1) After a careful examination of the simple model 
it is found to fail to describe adequately the pressure 
along the bottom and it is therefore not useful for future 
efforts in describing the pressure radiated into the 
POertom. 

(2) On the other hand the method of images gives a 
good prediction of pressure amplitude and phase distribu- 
tions along the interface and therefore offers us hope in 
attempting to predict the properties of the beam of sound 


mathe bottom. 


ay 





APPENDIX A 


A computer program for the calculation of the pressure 
amplitude and phase distribution along the bottom by the 


method of images. 
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Pine.) BaeogWdine SesuGSelh MY NOTH PALE OF IMAGES Pre 

E Se noe Pe eemiWoshele FAIK UF IMAGES AKE PyZ 

Gunner SES FAN 

JCP Cte iOnnielo ky LMAGE, APEX AND BOTTOM Fy R 

ov Pangie pene Unm CUNVJENTENCE Fir 

RK [OGRE ereCrIgnN LOSS OF CNT OTH IMAGE as A 

N Vet ieenoe ON OTH PAL OF LMAGES a 

M Me tieesOa hOUM BOUNCE rok 

MM WCU eNO cina OF LOPTOM BOUNCES pyre 


aoc Oeil y THE NUTS INiGe THE KECEITVING FOINT Fr 
AND CHE BOTTOM 


Lik iv eooe None lhe hiniaGe TOSTHE KECEIVING FOINT Ban 
NORMALIZED BY DIVIDING BY XO 
(ite (Ghetto rnomanaole TO Tre BUTTOM Fy R 


Ore NMEPCAneE REFLECTION) 
Cre ee Ciietr tte Lr ORAZING ANGLE EXEEEGS CRITICAL ANGLE Fy R 


eget POsiiittewetN REPEECTION COERRITELENT Paz 
REPL eee eee COR PE LEN TGF INTERMED TATE EOUNCE FxZ 
hee omniamee i POM CORPRICILENT OF THE LAST BOUNCE Pee 
PISCE Tee ouuURC EE ULRECT IVIg y Fat 
ro Pmivine tem FOR CONVENIENCE mars 


—— met ee cr cre cre cree ec ee ee ee eee ee ee eee ee ee ee eee ee eee ee eee ee eee ee 
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(CFP ROGKAM) ) 


Qgonn 


BUBROUTINE PRES (UISTXsPAHFsFHAS) 


C 
C DEFINE CHARACTERG OF VAKRIARLEsFARAMETERS AND CONSTANTS 
C 
Pe Clie ne AL eS CA-tlyO0—Z2)s INTEGER CI=N) 
CUpetetwroer online bls hKEPILNOyKr Ze PRES »F 
COMMON C2LsKC2QL»BETAsWN FAI 2s ANGLO » X stl F 
C 


ert Chibi ze sCOnMPEEX FRESSURE P=0.04+J0,0 
P=DCMFLX(O,00070.0L0) 


G RESET COUNTER N=0 
Cc 
N=0 
10 CONTINUE 
E t 
(C CALCULATE ANGLN ANI GETERMINE IF ANGLE IS IN THE RANGE OR NOT 
C fF ANGLE.GIr.FAI GO TO THE NEXT STEP 
c 
ANGLN=2,OL0XINTCC(Nt1L. 3/2.) KBETA4(-1) KXNKXANGLO 
IF (ANGLN.GE.FAI) GO TO SO 
C 
C CALCULATE FARAMETER THETAO AND DIR 
GC 
Lo=0COS CANGLN)-OISTX 
THETAO=ATANOCISINCANGIN) »f0) 
UR=CSQRT CL. OLOFDISTX**k2—-2,.0D0XDISTX*IICOS (ANGLN) ) 
c 
C INITIALIZE THE TOTAL REFLECTION COEFFCIENT R=1.0+40.0 
c 
R=ftCMFILX(1.0010°0,0L0) 
C 
C OEFINE MM BY EACH FATH WHICH IS DETERMINED BY THE N TH IMAGE 
c 
MM= IDENT CANGLN/(2.QUOXLETA) ) 
L=MMN+1 
00 40 M=tLeLl 
By 
C CALCULATE FARAMETERS THETAM AND CHECK 
c 
THE TAM=THETAO-2.OLOXIIFILOAT(CM-1)*BETA 
CHECK=1,000-C21 «2D *€0COS (CTHETAM) KX2 
e 
C IDENTIFY THE CNCLIENT ANGLE WHICH IS LESS THAN THRE CRITICAL 
C ANGLE OR NOT 
ce 
IF (CHECK.GT.0.000) GO TO 20 
c 
C CALCULATE FARAMETER FSAL FOR THE REFLECTION COEFFCIENT OF EACH 
C BOUNCE, THERE ARE TWO WAYS WHICH DEFEND ON THE LDOENTIFICATION 
C OF A CHECK 
e 
PSALT=TICMPLX(O,000»-0SURL(-CHECK) /DSINCTHETAM) ) 
GO TO 30 
20 CONTINUE 
FSAI=UCMPLXCOSQRT (CHECK) /DIIS INC THETAM) 20.000) 
30 CONTINUE 
i 
C CALCULATE REFL 
GE 


REPE=CRCei=P SALI CRO eT SAL) 
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i ww) oe ew) 


GaGa (1G) 


(PQ) BUS 


fee sOUNnG ees tHe EAS i rho THE SOURCE: THEN KEFLNO=REFL 
Pe CMistas i) hKEPLMO=REPL 
eeu towne TON REPEECTION €S8EFFECIENT R FOR EACH FATH 


R=RKXREFL 
1) CONTINUE 


POPeUeAte PAROMETEAS Z AND PRES 

ino ieto THe MiREET IYI Ty OF [He SOURCE 
eo evo ONLY WHEN THE SOURCE IS G@IRECT LONAL 
MITER WISE DinkET=140 


ZTDCMPLX(O0.000r-WNKXXKDR ) 
ULIRET=USINCO,OLODOXWNKDISIN( THE TAM-ANGLO))/ 
C CO.OLOUDO KWN¥DSING THE TAM-ANGLO) ) 
Pine os ine AGWExi (eo ke =D) KIMI CONT 1.) 72.) x 
te (1.000tL. OBO/REFLNO) «R/DR 


CALCULATE THE FRESSURE F WHICH [S A SUM OF FRES 
FSP+PRES 
SET TO THE NEXT FAIR OF TMAGE 


N=N+1 
CO. TO 10 
Ss0n CU ENUE 


Peewee ote mite ne TURN BACK TO THE MAIN ROUTINE 


Partials << Lei CF) 

PRE AL=F 

Pee = Real 
PP=PKICMELXCO,0T0,1.000> 
FIMAG=OUSIGN(CCUARS CE)» -FF) 
PRHAS=DAtTANetr LRAG rr REAL) 
RETURN 

ENTE 
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AM Sed NBME CAs: 


INCIDENT ANGLE FOR A DISTANCE SOURCE 


General expression of an incident angle is given 


by Eq. (II 39-a) 


6 = 9 > Pains’ Pisa) 


where the maximum number of bounces M is limited by 


M = INT al fr1-40) 
For a distance source 
Q ie (II-48) 


A few incident angles se for a distance source are shown 
in Table B-l. From this table it is seen that the reflec- 


men coefficient for the mee image is 


and by mathematical induction Eqs. II-50 can be deduced. 
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